Columnar and lamellar phases in attractive colloidal systems 
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In colloidal suspensions, at low volume fraction and temperature, dynamical arrest occurs via 
the growth of elongated structures, that aggregate to form a connected network at gelation. Here 
we show that, in the region of parameter space where gelation occurs, the stable thermodynamical 
phase is a crystalline columnar one. Near and above the gelation threshold, the disordered spanning 
network slowly evolves and finally orders to form the crystalline structure. At higher volume fractions 
the stable phase is a lamellar one, that seems to have a still longer ordering time. 
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In colloidal suspensions solid (or liquid) mesoscopic 
particles are dispersed in another substance. These sys- 
tems, like blood, proteins in water, milk, black ink or 
paints, are important in our everyday lives, in biology 
and industry [1, 2]. It is crucial, for example, to control 
the process of aggregation in paint and paper industries 
[3] , or to favour the protein crystallization in the produc- 
tion of pharmaceuticals and photonic crystals [4, 5]. 

A practical and exciting feature of colloidal suspen- 
sions is that the interaction energy between particles can 
be well controlled [6-8]. In fact particles can be coated 
and stabilized leading to a hard sphere behaviour, and 
an attractive depletion interaction can be brought out 
by adding some non-adsorbing polymers. The range and 
strength of the potential are controlled respectively by 
the size and concentration of the polymer [8, 9]. Re- 
cent experimental works highlighted the presence of a net 
charge on colloidal particles [7, 10] giving rise to a long 
range electrostatic repulsion in addition to the depletion 
attraction. 

The competition between attractive and repulsive in- 
teractions produces a rich phenomenology and a complex 
behavior as far as structural and dynamical properties 
are concerned. For particular choices of the interaction 
parameters, the aggregation of particles is favoured but 
the liquid-gas phase transition can be avoided and the 
cluster size can be stabilized at an optimum value [11]. 
Experimentally, such a cluster phase made of small equi- 
librium monodisperse clusters is observed using confocal 
microscopy at low volume fraction and low temperature 
(or high attraction strength) [7, 10, 12]. Increasing the 
volume fraction, the system is transformed from an er- 
godic cluster liquid into a nonergodic gel [10, 12], where 
structural arrest occurs. Using molecular dynamic sim- 
ulations, we showed that such structural arrest is cru- 
cially related to the formation of a long living span- 
ning cluster, providing evidence for the percolation na- 



ture of the colloidal gel transition at low volume fraction 
and low temperature [13, 14]. This scenario was con- 
firmed by recent experiments [10] and molecular dynam- 
ics simulations [15], where it was shown that increasing 
the volume fraction clusters coalesce into elongated struc- 
tures eventually forming a disordered spanning network. 
A realistic framework for the modelization of these sys- 
tems is represented by DLVO interaction potentials [16], 
which combine short-range attractions and long-range re- 
pulsions. A suited choice of DLVO models has in fact 
allowed to reproduce, by means of molecular dynamics 
calculations, many experimental observations, like the 
cluster phase and the gel-like slow dynamics [13-15]. 

On the other hand, competing interactions have been 
studied in many other systems, ranging from spin sys- 
tems to aqueous surfactants or mixtures of block copoly- 
mers, and often lead to pattern formation or to the cre- 
ation of periodic phases [17-23]. 

In this paper, we simulate by molecular dynamics a 
system composed of monodisperse particles, interacting 
with a short range attraction and a long range repul- 
sion in analogy with DLVO models, for a large range 
of temperatures and volume fractions. At low tempera- 
ture, increasing the volume fraction in the region of phase 
space where the system forms a percolating network and 
waiting long enough, we observe that the system spon- 
taneously orders, to form a periodic structure composed 
by parallel columns of particles. This finding strongly 
suggests that the transition to the gel phase, observed 
in experiments and in numerical simulations, happens 
in a "supercooled" region, namely in a disordered phase 
which is metastable with respect to crystallization. This 
is supported also by the results on a mean field model 
[21]. We then study the phase diagram of the system, by 
evaluating the free energy of the disordered and ordered 
phases, and find the region where the columnar phase is 
stable. We also locate the region of phase space where 



the stable phase is the lamellar one. Such phase does not 
form spontaneously within the observation times, there- 
fore indicating a much longer nucleation time. 

We have considered a system made of 1300 < N < 
2300 particles. The particles interact through the effec- 
tive interaction potential [14]: 
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where ai — 2.3, a,2 — 6, 03 = 3.5, and A = 2.5. The 
potential is truncated and shifted to zero at a distance 
of 3.5cr. The temperature T is in units of e/ks, where 
fcs is the Boltzmann constant. The number of particles 
varies as we consider different volume fractions, defined 
using the equivalent system of spheres of diameter a as 
4> = 4/3ira 3 N/L 3 in a simulation box of size L. We 
have performed Newtonian molecular dynamics at con- 
stant NVT using the velocity Verlet algorithm and the 
Nose-Hoover thermostat [24] with time step At = O.Olio 
(where to — y/ ma 2 /e and m is the mass of the particles). 
For a given volume fraction 0, the system is first thermal- 
ized at very high temperature and then quenched to the 
desired temperature. We let it evolve for times ranging 
between 10 5 t and 10 6 t (10 7 -10 8 MD steps). 

At low volume fraction (</> < 0.2) and high temper- 
ature (T > 0.25), the system remains in a disordered 
configuration during the simulation time window. De- 
creasing the temperature, we observe the formation of 
long living clusters with a typical size of about ten par- 
ticles. At temperature T = 0.25 and for volume frac- 
tions <f> between 0.15 and 0.23, we find that, shortly after 
the quench, these clusters aggregate into locally rod-like 
structures with a length from a few to about ten molec- 
ular diameters which finally coalesce into a disordered 
network. At longer times, between 5 x 10 4 to an d 10 6 to, 
these elongated rod-like structures start to spontaneously 
order in a two-dimensional hexagonal packing of parallel 
columnar structures, with a fast drop in the energy, as 
illustrated in Fig. 1(a). We also observe that locally par- 
ticles rearrange in a spiral structure which optimizes the 
competing interactions, very similar to the Bernal spiral 
[25], where each particle has six neighbors. Such a local 
structure is also indicated by experimental observations 
[10] and found in numerical simulations with similar in- 
teractions [15]. 

In order to infer the phase diagram of the system at low 
volume fraction, we first studied the state of the system 
at very low temperature. In that limit, entropy can be 
neglected and the equilibrium state of the system should 
be the crystalline structure with the lowest potential en- 
ergy. Of course, an exact determination of such a struc- 
ture, that may depend on the volume fraction, is far from 
trivial. We have instead selected a few possible structures 
and made a comparative analysis by calculating numer- 
ically their potential energy as a function of the volume 
fraction. The structures considered are the following (in 
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FIG. 1: a) Snapshot of one of the systems after the ordering 
(simulation time t = 10 J to), taken from a direction parallel to 
the columns, showing the two-dimensional arrangement of the 
axes. Since columns are not perfectly straight, their diameters 
appear to be larger than in reality; b) Lamellar phase at <j> = 
0.26 and T = 0.2. 



each case, the lattice spacing between elements is chosen 
to provide the desired volume fraction): 

1) Cluster crystal: Three-dimensional hexagonal close 
packing of nearly spherical clusters, each cluster being 
composed of thirteen particles (a central one in contact 
with twelve external ones). In this case particles have 
5.54 neighbors on the average. These clusters were cho- 
sen because of their similarity to the long-living clusters 
observed at low temperature and volume fraction, com- 
posed generally by one particle surrounded by eight to 
eleven neighbors. 

2) Columnar phase: Two-dimensional hexagonal packing 
of parallel Bernal spirals. In this case, each particle has 
six neighbors. 

3) Lamellar phase: Parallel planes, each plane being 
formed by one or more layers of hexagonally packed par- 
ticles (Fig. 1(b)). In this case, the average number of 
neighbors per particle is between six (one layer) and 
twelve (when the number of layers increases the structure 
becomes an hexagonal close-packed lattice of particles). 

After having placed the particles to form these struc- 
tures with the distance between nearest neighbors cor- 
responding to the minimum of the pair-potential, we 
have calculated the nearest minimum of the potential en- 
ergy landscape, by minimizing the total potential energy 
U(r±, . . . , r^r) as a function of the particle positions, us- 
ing a conjugate-gradient algorithm [27]. The potential 
energy per particle Uo/N of the minimum is shown in 
Fig.2. 

The cluster crystal never provides the lowest energy 
because of the limited number of neighbors per particle. 
For volume fractions tj> < 0.235, the columnar phase of 
Bernal spirals is energetically preferred. Indeed, in this 
case the combination of a relatively high number of first 
neighbors and a low number of second and third neigh- 
bors minimizes the repulsive part of the interaction po- 
tential. Increasing the volume fraction, the columns get 
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FIG. 2: Minimum of the total potential energy per particle, 
Uo /N, for several structures as a function of the volume frac- 
tion. 

closer and the repulsion between them becomes relevant. 
At <j> > 0.235, the potential energy of the columnar phase 
becomes higher than that of the two-layer lamellae. In 
fact, since lamellae are more distant than columns for a 
given volume fraction, the repulsion between lamellae is 
weaker. It is interesting to note that one-layer lamellae 
are never energetically preferred and that their minimum 
potential energy becomes close to the one of the columns 
at high volume fraction. At (j> > 0.42, the three-layer 
lamellae exhibit the lowest potential energy, and so on. 
Similar results were recently found in a 2d DLVO system 
[23]. 

In order to test the stability of these structures as the 
temperature increases, we performe the following exper- 
iment. Given a volume fraction and a temperature, we 
let evolve a crystalline configuration of columns or two- 
layer lamellae by means of molecular dynamics at con- 
stant temperature for 10 6 £o- The obtained phase diagram 
is shown in Fig. 3. Circles represent state points where 
both columns and lamellae break down before the end 
of the simulation (disordered phase). Triangles are the 
points where lamellae break down, but columns are stable 
(columnar phase). Squares are the points where lamel- 
lae are stable and columns break down (lamellar phase). 
Finally, the points where triangles and squares overlap 
correspond to the states where both columns and lamel- 
lae remain stable until t = 10 6 io- We note that near the 
boundaries with the disordered phase, columns or lamel- 
lae become fuzzy, but the modulation of the density is 
still clearly visible. 

To check the results found by molecular dynamics sim- 
ulations and to better estimate the phase boundaries, we 
have computed the free energy of the three phases that 



FIG. 3: T-cf> phase diagram of the system. Circles repre- 
sent the disordered phase, triangles the columnar phase and 
squares the lamellar phase. Stars represent points where the 
calculated free energy of the two phases crosses. Solid lines 
are a guide for the eyes. 

seem to be relevant for <p < 0.4: disordered phase, colum- 
nar phase, and two-layer lamellar phase. The free energy 
of the disordered phase is computed by thermodynamic 
integration [26] along the following path: from the per- 
fect gas limit (V — > oo) down to the desired volume V 
along the isotherm To = 1. The free energy at (V,To) is 
therefore given by 

oo 

NT \VJ 2 \mT J J NT 

v 

(2) 

where AP(V) — P(V) — -^P 1 is the excess pressure with 
respect to the perfect gas. Since values of h and the mass 
m merely shift the free energy by a constant, we set h = 1 
and m = 1. We have then integrated along the isochore 
V, from To to the desired temperature, and obtained 

^ = ^ + hJ^\J^ d T>, (3) 
NT NT 2 \T J J NT' 2 

T 

Starting from T = 0, we have calculated the free energy 
of the two crystalline phases by integrating along an iso- 
chore. At very low temperatures the potential energy 
landscape can be approximated by a parabolic function 

^ 3N 

U = U + -J2^l (4) 

i—l 

where are the eigenvalues of the Hessian matrix at 
the minimum of the potential. The free energy is then 
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FIG. 4: Free energy per particle F(T)/N of the disordered 
phase (circles) and of the columnar phase (triangles) as a 
function of temperature at tf> = 0.1308. Inset: potential en- 
ergy per particle U (T) /N. 



calculated as 
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where AU(T) = U(T) - l/ - MZ. 

In Fig. 4 the free energy of the disordered phase and of 
the columnar phase at <f> — 0.1308 are shown as a function 
of the temperature. The curves cross at the temperature 
(T c = 0.17) corresponding to the first order transition 
from the columnar phase to the disordered phase. This, 
together with the other transition points, are marked in 
Fig. 3 by a star [28]. 

These results suggest the following scenario for col- 
loidal gelation at low volume fraction and low temper- 
ature. Just remind that in real colloidal systems the 
control parameter corresponding to low temperature is 
a high effective attraction strength. The competition 
between attractive and repulsive interaction induces a 
typical modulation length clearly detected at very low 
volume fraction in the cluster phase. As the volume frac- 
tion is increased, clusters are prone to aggregate in spite 
of the long-range repulsion. This results in the grow- 
ing of elongated structures which keep track of the mod- 
ulation length. Eventually, a modulated phase forms, 
with a columnar geometry or a lamellar one at higher 
volume fraction. Of course the precursors of the modu- 
lated structures will have defects such as local inhomo- 
geneities and branching points. Furthermore, the slow 
dynamics due to the viscosity of the solvent and to the 
imperfect shape of the particles may also hinder the or- 
dered phase, producing long-living metastable disordered 



states. Hence, when the volume fraction is high enough, 
an interconnected spanning structure is formed and a gel- 
like behavior is observed as it is frequently the case in the 
experiments. 
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